Erratum to: Molecular modelling study of 2-phenylethynyladenosine (PEAdo) derivatives as highly selective A3 adenosine receptor ligands

A series of 2-phenylethynyladenosine (PEAdo) derivatives substituted in the N6- and 4′position was synthesised and the new derivatives were tested at the four human adenosine receptors stably transfected into Chinese hamster ovary (CHO) cells, using radioligand binding studies (A1, A2A, A3) or adenylyl cyclase activity assay (A2B). Binding studies showed that the presence of a phenyl ethynyl group in the 2 position of adenosine favoured the interaction with A3 receptors, resulting in compounds endowed with high affinity and selectivity for the A3 subtype. Additional substitution of the N6- and 4′position increases both A3 affinity and selectivity. The results showed that the new compounds have a good affinity for the A3 receptor and in particular, the N6-methoxy-2-phenylethynyl-5′N-methylcarboxamidoadenosine, with a Ki at A3 of 1.9 nM and a selectivity A1/A3 and A2A/A3 of 4,800- and 8,600-fold, respectively. Therefore, it is one of the most potent and selective agonists at the human A3 adenosine receptor subtype reported so far. Furthermore, functional assays of inhibition of 10 μM forskolin-stimulated cAMP production via the adenosine A3 receptor revealed that the new trisubstituted adenosine derivatives behave as full agonist of this receptor subtype. Docking analysis of these compounds was performed at a homology model of the human A3 receptor based on the bovine rhodopsin crystal structure as template, and the results are in accordance with the biological data

PubChem3D: a new resource for scientists

<p>Abstract</p> <p>Background</p> <p>PubChem is an open repository for small molecules and their experimental biological activity. PubChem integrates and provides search, retrieval, visualization, analysis, and programmatic access tools in an effort to maximize the utility of contributed information. There are many diverse chemical structures with similar biological efficacies against targets available in PubChem that are difficult to interrelate using traditional 2-D similarity methods. A new layer called PubChem3D is added to PubChem to assist in this analysis.</p> <p>Description</p> <p>PubChem generates a 3-D conformer model description for 92.3% of all records in the PubChem Compound database (when considering the parent compound of salts). Each of these conformer models is sampled to remove redundancy, guaranteeing a minimum (non-hydrogen atom pair-wise) RMSD between conformers. A diverse conformer ordering gives a maximal description of the conformational diversity of a molecule when only a subset of available conformers is used. A pre-computed search per compound record gives immediate access to a set of 3-D similar compounds (called "Similar Conformers") in PubChem and their respective superpositions. Systematic augmentation of PubChem resources to include a 3-D layer provides users with new capabilities to search, subset, visualize, analyze, and download data.</p> <p>A series of retrospective studies help to demonstrate important connections between chemical structures and their biological function that are not obvious using 2-D similarity but are readily apparent by 3-D similarity.</p> <p>Conclusions</p> <p>The addition of PubChem3D to the existing contents of PubChem is a considerable achievement, given the scope, scale, and the fact that the resource is publicly accessible and free. With the ability to uncover latent structure-activity relationships of chemical structures, while complementing 2-D similarity analysis approaches, PubChem3D represents a new resource for scientists to exploit when exploring the biological annotations in PubChem.</p

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation

Average and extreme multi-atom Van der Waals interactions: Strong coupling of multi-atom Van der Waals interactions with covalent bonding

<p>Abstract</p> <p>Background</p> <p>The prediction of ligand binding or protein structure requires very accurate force field potentials – even small errors in force field potentials can make a 'wrong' structure (from the billions possible) more stable than the single, 'correct' one. However, despite huge efforts to optimize them, currently-used all-atom force fields are still not able, in a vast majority of cases, even to keep a protein molecule in its native conformation in the course of molecular dynamics simulations or to bring an approximate, homology-based model of protein structure closer to its native conformation.</p> <p>Results</p> <p>A strict analysis shows that a specific coupling of multi-atom Van der Waals interactions with covalent bonding can, in extreme cases, increase (or decrease) the interaction energy by about 20–40% at certain angles between the direction of interaction and the covalent bond. It is also shown that on average multi-body effects decrease the total Van der Waals energy in proportion to the square root of the electronic component of dielectric permittivity corresponding to dipole-dipole interactions at small distances, where Van der Waals interactions take place.</p> <p>Conclusion</p> <p>The study shows that currently-ignored multi-atom Van der Waals interactions can, in certain instances, lead to significant energy effects, comparable to those caused by the replacement of atoms (for instance, C by N) in conventional pairwise Van der Waals interactions.</p

Fast 3D shape screening of large chemical databases through alignment-recycling

<p>Abstract</p> <p>Background</p> <p>Large chemical databases require fast, efficient, and simple ways of looking for similar structures. Although such tasks are now fairly well resolved for graph-based similarity queries, they remain an issue for 3D approaches, particularly for those based on 3D shape overlays. Inspired by a recent technique developed to compare molecular shapes, we designed a hybrid methodology, alignment-recycling, that enables efficient retrieval and alignment of structures with similar 3D shapes.</p> <p>Results</p> <p>Using a dataset of more than one million PubChem compounds of limited size (< 28 heavy atoms) and flexibility (< 6 rotatable bonds), we obtained a set of a few thousand diverse structures covering entirely the 3D shape space of the conformers of the dataset. Transformation matrices gathered from the overlays between these diverse structures and the 3D conformer dataset allowed us to drastically (100-fold) reduce the CPU time required for shape overlay. The alignment-recycling heuristic produces results consistent with <it>de novo </it>alignment calculation, with better than 80% hit list overlap on average.</p> <p>Conclusion</p> <p>Overlay-based 3D methods are computationally demanding when searching large databases. Alignment-recycling reduces the CPU time to perform shape similarity searches by breaking the alignment problem into three steps: selection of diverse shapes to describe the database shape-space; overlay of the database conformers to the diverse shapes; and non-optimized overlay of query and database conformers using common reference shapes. The precomputation, required by the first two steps, is a significant cost of the method; however, once performed, querying is two orders of magnitude faster. Extensions and variations of this methodology, for example, to handle more flexible and larger small-molecules are discussed.</p

Confab - Systematic generation of diverse low-energy conformers

<p>Abstract</p> <p>Background</p> <p>Many computational chemistry analyses require the generation of conformers, either on-the-fly, or in advance. We present Confab, an open source command-line application for the systematic generation of low-energy conformers according to a diversity criterion.</p> <p>Results</p> <p>Confab generates conformations using the 'torsion driving approach' which involves iterating systematically through a set of allowed torsion angles for each rotatable bond. Energy is assessed using the MMFF94 forcefield. Diversity is measured using the heavy-atom root-mean-square deviation (RMSD) relative to conformers already stored. We investigated the recovery of crystal structures for a dataset of 1000 ligands from the Protein Data Bank with fewer than 1 million conformations. Confab can recover 97% of the molecules to within 1.5 Å at a diversity level of 1.5 Å and an energy cutoff of 50 kcal/mol.</p> <p>Conclusions</p> <p>Confab is available from <url>http://confab.googlecode.com</url>.</p

CheS-Mapper - Chemical Space Mapping and Visualization in 3D

Analyzing chemical datasets is a challenging task for scientific researchers in the field of chemoinformatics. It is important, yet difficult to understand the relationship between the structure of chemical compounds, their physico-chemical properties, and biological or toxic effects. To that respect, visualization tools can help to better comprehend the underlying correlations. Our recently developed 3D molecular viewer CheS-Mapper (Chemical Space Mapper) divides large datasets into clusters of similar compounds and consequently arranges them in 3D space, such that their spatial proximity reflects their similarity. The user can indirectly determine similarity, by selecting which features to employ in the process. The tool can use and calculate different kind of features, like structural fragments as well as quantitative chemical descriptors. These features can be highlighted within CheS-Mapper, which aids the chemist to better understand patterns and regularities and relate the observations to established scientific knowledge. As a final function, the tool can also be used to select and export specific subsets of a given dataset for further analysis

Molecular dynamics simulations and in silico peptide ligand screening of the Elk-1 ETS domain

Background: The Elk-1 transcription factor is a member of a group of proteins called ternary complex factors, which serve as a paradigm for gene regulation in response to extracellular signals. Its deregulation has been linked to multiple human diseases including the development of tumours. The work herein aims to inform the design of potential peptidomimetic compounds that can inhibit the formation of the Elk-1 dimer, which is key to Elk-1 stability. We have conducted molecular dynamics simulations of the Elk-1 ETS domain followed by virtual screening. Results: We show the ETS dimerisation site undergoes conformational reorganisation at the a1b1 loop. Through exhaustive screening of di- and tri-peptide libraries against a collection of ETS domain conformations representing the dynamics of the loop, we identified a series of potential binders for the Elk-1 dimer interface. The di-peptides showed no particular preference toward the binding site; however, the tri-peptides made specific interactions with residues: Glu17, Gln18 and Arg49 that are pivotal to the dimer interface. Conclusions: We have shown molecular dynamics simulations can be combined with virtual peptide screening to obtain an exhaustive docking protocol that incorporates dynamic fluctuations in a receptor. Based on our findings, we suggest experimental binding studies to be performed on the 12 SILE ranked tri-peptides as possible compounds for the design of inhibitors of Elk-1 dimerisation. It would also be reasonable to consider the score ranked tri-peptides as a comparative test to establish whether peptide size is a determinant factor of binding to the ETS domain